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1. Introduction 

The study of QSO clustering from large homogeneous surveys is likely to pro- 
vide unique information on the nature of large-scale structure (LSS) in the Uni- 
verse. Even with relatively modest aperture telescopes, QSOs are readily de- 
tectable from the present cosmic epoch to a time when the Universe was less 
than 25% of its present age. This approximately eight billion year time interval 
provides a unique baseline over which to study the evolution of structure in the 
Universe with cosmic time. QSO clustering also yields a direct measure of cluster- 
ing on comoving scales S> 100h -1 Mpdj], providing key information (see figure 1) 
on LSS at scales between those provided by the COBE results (> 1000h -1 Mpc, 
Bennett et al. 1994) and galaxy redshift surveys (< lOOh^Mpc, Colless et al, 
this volume) 

To date, QSO clustering has remained largely unexploited for LSS studies. This 
has largely been due to the lack of appropriate QSO catalogues. Although over 
10000 QSOs with measured redshifts are now known (Veron & Veron-Cetty 1997, 
see figure 2) less than 20 % of these QSOs form part of homogeneous catalogues 
suitable for detailed clustering analysis. The goal of the 2dF QSO redshift survey, 
currently underway at the Anglo- Australian Telescope (AAT) , is to provide a new 
homogeneous catalogue of almost 30000 QSOs with which to carry out such an 
analysis. The survey is almost two orders of magnitude larger than the previous 
largest homogeneous survey to this depth (Boyle et al. 1990). 

In § 2 we review the existing results on QSO clustering and we describe the 
creation of the 2dF QSO redshift survey input catalogue in § 3. We report on the 
initial results from the survey in § 4 and discuss the future prospects for this and 
other QSO surveys in § 5. 



2. QSO Clustering 



f Throughout this paper we use comoving distances. r comov i ng = (1 + z)r propcT , and define the Hubble 
constant as Hq = lOOhkm s — 1 Mpc _1 
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Figure 1. The QSO power spectrum, P(k), based on a Zel'dovich approximation in a 
(2250h _1 Mpc) 3 box containing 192 3 particles with z max = 2.2. A rCDM model with Q = 1, 
r = 0.2 as = 0.67 and QSO bias 6q = 1.8 at z = is assumed. Errors on P(k) are based on the 
rms scatter from five simulations. Note the 2dF QSO redshift survey comprises two such slices, 
effectively reducing errors by 1.4. Simulation taken from Croom et al. (1998, in preparation). 
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Figure 2. The all-sky distribution of the ~ 10000 known QSOs from the catalogue of Veron & 
Veron-Cetty (1997). The thirty UKST fields of the 2dF QSO redshift survey are also shown in 
their two declination strips. 
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Figure 3. The QSO auto-correlation function £c>(r) f° r the redshift range 0.3 < z < 2.2 derived 
by Croom & Shanks (1996) from the Durham/AAT (Boyle et al. 1990), CFHT (Crampton et 
al. 1989) and LBQS (Hewett et al. 1995) surveys. The total sample comprises ~ 1500 QSOs. 

(a) QSO- QSO clustering 

First attempts to detect QSO clustering from samples of QSOs yielded no sig- 
nificant detections of QSO clustering on any scale up to 1500h _1 Mpc (see e.g. 
Osmer 1981). Nevertheless, a number of studies (Osmer 1981, Webster 1982) 
did comment on the possibly genuine nature of specific clusters of QSOs at in- 
termediate scales ~ 50 — 100h _1 Mpc. The first positive detection of clustering 
amongst QSOs at small scales (< 10h _1 Mpc) was made by Shaver (1984) using 
an earlier version of the heterogeneous Veron & Veron-Cetty (1997) catalogue, 
which then comprised 2000 QSOs. This was followed by similar results based on 
more homogeneous QSO catalogues, comprising a few hundred QSOs (Shanks, 
Boyle & Peterson 1988, Iovino & Shaver 1988). The strength of QSO clustering in 
these studies was found to be similar to that of galaxies at the present day, with 
the two-point QSO auto-correlation function, £q(7"), fitted by a —1.8 power law 
£Q( r ) = ('"A'o) -1 ' 8 with a scale length ro ~ 6h _1 Mpc. None of these early QSO 
clustering studies found any significant QSO clustering on scales > 10h _1 Mpc. 

Suitable QSO sample sizes have only increased slowly over the past ten years, 
and latest analyses of QSO clustering (see e.g. Shanks & Croom 1996) still rely 
on a compilation of a few homogeneous catalogues (Durham/AAT: Boyle et al. 
1990, CFHT: Crampton et al. 1989, LBQS: Hewett et al. 1995) comprising less 
than 2000 QSOs in total. 

Using the Durham/AAT, CFHT and LBQS catalogues, Croom & Shanks (1996) 
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demonstrated that QSO auto-correlation function at scales less than 30h Mpc 
was best-fit by f Q (r) = ( r /r )- L59±0 - 46 , with r = 4.0 ± l.lh _1 Mpc (figure 3). 
Croom Sz Shanks (1996) also found no evidence for any evolution in the strength 
of clustering with redshift. By investigating the predications of a standard bias- 
ing model, they were able to demonstrate that the low values of vq obtained for 
z < 0.2 Seyfert galaxies (Boyle & Mo 1993) and IRAS-selected AGN (Georgan- 
topoulos & Shanks 1994) favoured either a low 0,q universe, CIq < 0.1, with low 
values for the galaxy bias parameter (b g ~ 1) or an = 1 universe with a higher 
bias parameter (6 g ~ 2). On the basis of this model, Croom & Shanks were also 
able to rule out a universe with f^o = 1, b g = 1 at the ~ 3a level. 

These results on the evolution of the QSO correlation function are, however, 
contested by other authors. In their earlier analysis, Iovino & Shaver (1988) found 
that the correlation length decreased towards high redshift, while more recently 
La Franca et al. (1997) claim to see an increase in the QSO correlation length 
with redshift, with r = 4.2h~ 1 Mpc for z < 1.4 QSOs and r = 9.1h _1 Mpc for 
z > 1.4 QSOs. 

The lack of agreement between these results serves to underline the poor statis- 
tics on which the detection of QSO clustering at small scales is currently based. 
There are less than 40 pairs of QSOs with comoving separations r < 10h~ 1 Mpc 
identified in such clustering studies, even when averaged over the full range in 
redshift 0.3 < z < 3. 

At intermediate scales, 30h _1 Mpc< r < 100h _1 Mpc, there have been a number 
of claims for marginal (2a) detections of features in the QSO auto-correlation 
function. However, none of these has stood up to critical examination with larger 
data-sets (see Croom & Shanks 1996). At these scales, the growth of LSS largely 
follows linear theory and thus features in £q(?") at these comoving separations 
should not change comoving scale as a function of redshift. If detected, such 
features would therefore provide a useful cosmological 'standard rod' with which 
to measure the deceleration parameter qo or the cosmological constant Ao- 

Croom and Shanks (1996) also compared £q(?*) i n t ne 8—50 h _1 Mpc region with 
an extrapolation of the CDM power spectrum normalised by the COBE data. This 
directly compares the fluctuations in the QSO distribution with those in the mass 
distribution and is thus a potentially powerful way to measure directly the QSO 
bias parameter. For a standard CDM model with f2o = 1, they obtained a QSO 
bias of 6q(z = 0) = 1.40^43 • For a CDM model with a cosmological constant 
(fi = 0.2, A = 0.8, r = 0.2 ~ fi h), they obtained b Q (z = 0) = 1.20±g }|. 

The inconsistency between the small scale results, which favour a high bias for 
J7o = 1) and the COBE comparison, which favour a low bias for £\) = 1, could 
be taken to mean that a low value for m ay be more appropriate. However, 
the errors are still large, and it will take surveys such as the 2dF QSO survey 
to reduce them to the point where more effective discrimination between these 
models is possible. 

Phillips et al. (1996) and Ballinger et al. (1997) have also developed a poten- 
tially useful geometric test which may be employed at these scales to place limits 
on the value of Ao- The test relies of the fact that, on average, QSO clusters 
should be spherically symmetric i.e. with as large an extent along the line of 
sight as transverse to it. Since Ao has a different effect on these two distances, 
clusters will appear to be spherically symmetric only when the correct value of Aq 
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is used. Ballinger et al. (1996) demonstrated that the 2dF QSO redshift survey 
would indeed be able to measure this effect, although effects such as redshift-space 
distortions serve to weaken the signal and complicate the analysis. 

At the largest scales, all measures of the correlation function give a signal con- 
sistent with zero clustering. At comoving separations greater than 100h _1 Mpc 
Croom & Shanks (1996) find £ Q (r) = ± 0.025. Isolated examples of QSO 'su- 
perclusters' with dimensions ~ 100 — 200h~ 1 Mpc have previously been identified 
(Crampton et al. 1991, Clowes & Campusano 1993), but it is unclear how repre- 
sentative they are of the QSO population as a whole. 

(b) QSO-galaxy clustering 

A limiting factor in the study of LSS with QSOs is the unknown bias factor 
which exists between the QSO population and the underlying mass distribution. 
Studies of the galaxy clustering around QSOs provide a useful constraint on the 
relationship between the QSO and galaxy bias. 

The study of the galaxy environment around QSOs was first put on a sound 
statistical footing by Stockton (1982) and later developed by Yee and his co- 
workers in a series of papers, see Ellingson et al. (1991) and references therein. 

The work of Ellingson et al. (1991) largely focussed on radio-loud QSOs, and 
demonstrated that their galaxy environment undergoes a dramatic increase in 
richness from low redshifts to intermediate redshifts, corresponding to an evo- 
lution in the galaxy clustering scale length around QSOs from ro ~ 5h~ 1 Mpc 
(consistent with field galaxies) at z ~ 0.3 to ro ~ 15h _1 Mpc (similar to Abell 
richness class 1) at z ~ 0.6. Evidence for similar levels of strong galaxy clustering 
around radio-loud QSOs at yet higher redshifts [z ~ 1) have been found by Tyson 
(1986) and Hintzen et al. (1991). 

In contrast, the environments of the vast majority of QSOs which are radio- 
quiet appear to significantly poorer. Smith et al. (1996) demonstrated that the 
QSO-galaxy cross-correlation function for a sample of over 100 X-ray selected 
radio-quiet AGN at z < 0.3 from the Einstein Medium Sensitivity Survey was 
identical to the galaxy auto-correlation function over a similar range in redshift. 

Smith (1998) has now extended this analysis to higher redshift QSOs in the 
same X-ray selected sample, finding evidence for the evolution in the QSO-galaxy 
cross-correlation function in the opposite sense to that witnessed for radio-loud 
QSOs. The QSO-galaxy scale length decreases from ro = 5h~ 1 Mpc at z < 0.3 
to ro ~ 2h~ 1 Mpc in the redshift range 0.5 < z < 0.7 (figure 4). The inferred 
evolution of the QSO-galaxy cross-correlation function may be expressed in the 
form £qg(0 oc r~ L8 (l+z)~ 3+e , with e ~ 0.8, equivalent to the growth of structure 
predicted by linear theory in an Qq = 1 Universe (Peacock 1997). 

This evolution is also similar to the evolution of the galaxy scale length as 
inferred indirectly from the galaxy angular auto-correlation function derived from 
deep imaging surveys (see Efstathiou 1995) and measured directly from deep 
redshift surveys (Le Fevre et al. 1996). 

Based on these results we can infer that radio-loud and radio-quiet QSOs in- 
habit different galaxy environments, probably indicative of different levels of bias 
between the two populations. This is further borne out by the QSO correlation 
length (ro ~ 18h _1 Mpc) inferred from two-dimensional radio-source catalogues 
(Loan et al. 1997). This level of clustering is more consistent with recent esti- 
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Figure 4. The QSO-galaxy cross- correlation functions, «j qg (#), for samples of X-ray selected 
QSOs at z < 0.3 (top), 0.3 < z < 0.5 (middle), and 0.5 < z < 0.7 (bottom). In the z < 0.3 
sample the w qg (#) is directly compared with the galaxy auto-correlation function (dashed line). 
In the higher redshift samples predicted cross-correlations functions are shown for QSO-galaxy 
scale lengths of r = 2h~ 1 Mpc and ro = 5.5h~ 1 Mpc under certain assumptions regarding the 
evolution of clustering. See Smith et al. (1996) and Smith (1998). 
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Figure 5. The 2dF QSO survey input catalogue based on the selection of stellar objects from 
UBjR UK Schmidt photographic material. Holes in the data-set represent areas excised due to 
bright stars 



mates of the galaxy cluster correlation length (ro ~ 14h _1 Mpc, Dalton et al. 
1994). However, the extent to which the radio-quiet QSO population is repre- 
sentative of the galaxy population will probably only be answered by detailed 
redshift-dependent QSOs clustering studies in the 2dF QSO redshift survey. 



3. The 2dF input catalogue 

(a) Creation 

The identification of the QSO candidates for the 2dF redshift survey was based 
on broadband U BjR colour selection from APM measurements of UK Schmidt 
(UKST) photographic plates. The survey area comprises 30 UKST fields, ar- 
ranged in two 75° x 5° declination strips centred on 5 = —30° and 5 = 0° (see 
figure 5). The 5 = —30° strip extends from a = 21 h 40 to a = 3 h 15 and the 
equatorial strip from a = 9 h 50 to a = 14 h 50. The total survey area is 740 deg 2 , 
when allowance is made for regions of sky excised around bright stars, and is a 
subset of the survey area for the 2dF galaxy redshift survey. 

In each UKST field APM measurements of one J plate, one R plate and up to 
four U plates/films were used to generate a catalogue of stellar objects with Bj < 
20.85. A sophisticated procedure was devised to ensure catalogue homogeneity 
across plate boundaries (Smith 1998). The criteria for inclusion in the catalogue 
were (U-Bj) < -0.36; (U-Bj) < 0.12 - 0.8(5,/ - i?); (Bj-R) < 0.05. The full 
C/j5ji?-selected catalogue is shown in figure 5. Based on the colours of QSOs which 
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Figure 6. The auto-correlation function for the 2dF QSO survey input catalogue. The curves 
denote the expected angular correlation function assuming ro = 6h Mpc (dashed line) or 
ro = 18h _1 Mpc (dot-dashed) for a —1.8 power law spatial correlation function and stable 
clustering in comoving space. 



have previously been identified in the survey region, we estimate the catalogue 
is ~ 93% complete for z < 2.2 QSOs and comprises over 55% QSOs (see Croom 
1997). Preliminary indications from the 2dF survey confirm this QSO fraction, 
with the principal contamination arising from galactic subdwarfs and compact 
blue galaxies. 

(b) Results 

Smith (1998) has measured the angular auto-correlation function, wqq(9) from 
the 2dF QSO input catalogue. To minimise the contamination from sources other 
than QSOs, Smith (1998) restricted the auto-correlation analysis to objects in the 
input catalogue with U — Bj < —0.5. Based on the derived contamination rates, 
this restricted sample comprises almost 80% QSOs. 

The angular auto-correlation function derived by Smith (1998) is plotted in 
figure 6. Smith (1998) reports a 3a detection of clustering consistent with a scale 
length of ro = 4.9±1.0h~ 1 Mpc, assuming the scale length of the QSO correlation 
remains constant in comoving space as a function of redshift. This is the first 
detection of QSO clustering on the plane of the sky from an optically-selected 
catalogue and the amplitude is consistent with the scale length of QSO clustering 
previously derived from spatial clustering studies. It is inconsistent with the much 
greater scale length of ro ~ 18h _1 Mpc inferred for radio-loud QSOs by Yee & 
Green (1987). 
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Figure 7. The number-redshift relation for the QSOs identified to date (March f998) in the 

2dF QSO redshift survey. 

At angular scales which are multiples of the mean field-to-field separation, 5.5° 
and 11°, wqq(0) does show systematic deviations from zero amplitude. Although 
these deviations are small (< 0.005) they may represent the fundamental limit 
of the accuracy to which large-scale structure can be measured at scales r ~ 
400h~ 1 Mpc (the comoving separation corresponding to 5° at z = 1.5). 



4. The 2dF QSO redshift survey 

2dF spectroscopic observations for the QSO redshift survey began in January 
1997, although the bulk of the redshifts have been obtained in the October 1997 
and January 1998 observing runs as the 2dF system has gained functionality 
(increased number of fibres, faster field re-configuration times). Over 1000 QSO 
redshifts have been obtained with the 2dF, and the 2dF QSO redshift survey is 
already the largest single homogeneous QSO catalogue in existence. Initial results 
from this survey are shown in figures 7-10. 

Figure 7 shows the number-redshift, n(z), relation for the sample. The n(z) 
shows a broad flat distribution, with no obvious redshift peaks or troughs due 
to selection biases. The auto-correlation function plotted in figure 8 is consistent 
with the form previously derived by Croom & Shanks (1996); £(r) = (r/4h -1 )~ L6 . 
Although the error bars are still large, there is still no strong evidence for evo- 
lution of the correlation length with redshift (figure 9). The predictions for two 
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Figure 8. The spatial auto-correlation function for QSO identified to date in the 2dF QSO 
redshift survey. The dotted line denotes £q(t) = (r/4h _1 Mpc) -1 ' 6 

linear theory models of structure growth confirm the result of Croom & Shanks 
(1996) that the lack of evolution favours a low O universe in such models. 

With 97% of the 2dF QSO redshift survey remaining to be completed, the 
survey cone diagram plotted in figure 10 is currently rather empty, but it does 
serve to illustrate the eventual scale of the 2dF QSO redshift survey. 

The input catalogue is also useful as a resource for other QSO studies. In a sep- 
arate project, we have cross-correlated the input catalogue with the NRAO/VLA 
sky survey (NVSS, Condon et al. 1998), yielding over 600 detections (correspond- 
ing to 5% of the QSOs in the input catalogue). Keck spectra have been obtained 
for almost 100 of these sources, the vast majority of which (> 95%) are QSOs. 

We have also begun a more conventional long-slit spectroscopic program to 
obtain redshifts for all input catalogue pairs with angular separations less than 
20arcsec. Such objects are candidate wide-separation gravitational lenses or bi- 
nary QSOs. The incidence of wide-separation gravitational lenses can be used 
to place constraints on the value of the cosmological constant. Following Maoz 
et al. (1997), Croom (1997) calculates that lensing by 'typical' cluster masses 
(2 x 10 14 M©) would result in between four (S7o = 1, Ao = 0) and fourteen 
(S7o = 0.4, Ao = 0.6) 5 — 20-arcsec QSO lenses being seen in the survey. To date, 
we have surveyed approximately 25% of the catalogue pairs at these separations, 
and found only one possible lens candidate (figure 11). The two components of 
the possible lens (Q0101.8-3012A,B) have a separation of 17arcsec on the sky. 
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Figure 9. £q(10), the mean amplitude of the QSO correlation function below 10h _1 Mpc, plotted 
as a function of redshift (filled squares). Estimates are based on the combined Durham/ A AT, 
CFHT, LBQS and the 2dF survey data currently available. The point at z ~ (open square) 
is derived from the AGN correlation functions of Boyle & Mo (1993) and Georgantopoulos & 
Shanks (1994). Model predictions for the growth of structure based on linear theory in an Q = 1 
(solid line) and fi = 0.1 Universe (dashed line) are shown. 

At the redshift of the system (z = 0.89), this corresponds to a projected proper 
spatial separation of 71h~ 1 kpc. 



5. Future Prospects 

Currently the 2dF QSO redshift survey is only 3% complete. Yet it is already 
the single largest homogeneous QSO redshift survey yet compiled. Many of the 
cosmological tests proposed above will become possible when the survey is com- 
pleted, providing measures of Qq, 6q and Ao as well as new information on the 
redshift evolution of structure in the Universe and on LSS at scales of many 
hundreds to thousands of megaparsecs. In addition the catalogue will provide a 
unique resource for damped Lya systems, and common absorption line systems. 
It may also prove possible to derive information on the power spectrum of mass 
fluctuations at z > 2 through the Lya forest absorption via the method of Croft 
et al. (1997). 

The next decade will see the completion of the 2dF QSO survey together with 
the Sloan Digital Sky Survey (10 5 QSOs; see Margon, this volume). The LAMOST 
survey (10 6 QSOs, see Chu & Zhao 1997) is also planned to be well into its 
operational phase. When QSOs were first discovered, the hope was they they 
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Figure 10. Cone plots for QSOs identified in the 0° and —30° declination strips from 2dF 

observations to date. 
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Figure 11. Spectra for the QSO pair Q0101.8-3012A,B, with the ratio spectrum plotted at the 

bottom. 

could be used to carry out standard cosmological tests. Now, almost forty years 
on, with the advent of these new surveys, this goal may finally be realised. 
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